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ABSTRACT 



DUDLEY KNOX DBBARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA !?3940 



The design of second order sampled analog filters using 
Z-transform techniques developed in digital theory was studied. 
A second order recursive CTD- filter was implemented and 
investigated. Deviation from theoretical frequency 
response was found to be partly due to frequency dependent 
circuitry. 

MTI- simulation was performed, and the ability of the 
CTD-filter to cancel clutter and pass doppler frequencies 
was demonstrated. Unexplained glitches in the filter 
output were noted. 
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I. INTRODUCTION 



The field of Signal Processing is expanding rapidly, 
indicating an increasing importance. Today, digital signal 
processing is being extensively explored and employed. 

High speed digital computers with necessary peripheral 
equipment such as analog-to-digital and digital- to-analog 
converters are required. Traditional analog systems as radar 
and sonar are going digital. There exists a wealth of other 
digital signal processing applications. 

For some applications, however, the use of digital 
processing may be termed "overkill", even if there is a 
distinct requirement for the inherent advantages of digital 
processing over analog processing. 

A fairly new device technology, the Charge Transport 
Device technology (CTD's) offers new ways of signal processing 
which could be the answer for such cases. 

Charge Transport Devices refer to several families of devices; 
Bucket Brigade Devices (BBD's), Charge Coupled Devices (CCD’s), 
and Reticon Devices (i.e. SAD-100 and TAD-12) which are 
described in reference {3}. 

The unique feature of CTD's is that the signal is 
represented by discrete charge packets, which are moved 
through the device at a speed determ.ined by the clock 
frequency. The delay of the device is thus controlled by the 
clock frequency, and may be easily and momentarily changed. 

An excellent review of how the CCD works is given in 
reference {8}. 

CTD's find wide applications in imaging (both low light, 
normal light and infrared), analog time delay, multiplexing, 
transversal filtering, with fixed and variable v/eighting, 
recursive filtering, correlators, chirp-z transformers and 
so on. The general advantages of CTD signal processors 
over conventional digital processing are lov/er cost, power, 
weight and size. 



12 



Comb filters are described by their periodic transfer 
characteristics in the frequency domain. They can be classi- 
fied into two general types. The first type is the bandstcp 
or canceller type shown in Table II. It has strong attenu- 
ation in a narrow neighbourhood of a series of periodically 
separated frequencies and good transmission in between. 

The second type is the bandpass, or integrator type (Table II). 
It has good transmission in a narrow neighbourhood of a 
series of periodically separated frequencies and strong 
attenuation in between. All possible implementation 
techniques for such comb filters use some kind of delay de- 
vice . 

In the analog case, a quartz delay line has been 
extensively used. It is, hov;ever, temperature dependent 
and results have in general not been very satisfactory. 

For the digital case the delay device is the shift register, 
in V7hich the amount of delay is determined by the employed 
clock frequency. 

Today a new family of comb filters is being developed, using 

the new CTD's. In these, analog signals are first sampled 

and then delayed. As a result of the sampling, these devices 

not only have standard analog properties but also some 

digital properties such as aliasing and stability of delay. 

A CTD consists of N serial delay stages. A signal V7ill then 

be delayed N clock periods before it is fed back to the 

input in the recur’sive case. The frequency of recursive 

operation, f , is then only where f is the clock 

c c 

frequency, or sampling frequency. Therefore, in radar terms, 

f^ is equal to the Pulse Repetition Frequency (PRF), and 

there are N teeth of the comb within the Nyquist frequency 
2 

range from 0 to f^^/2. 

So far recursive filtering using CCD’s have not been 
extensively studied in the industry, since it was believed 
that it can only be done with extensive external electronic 
support such as fast and accurate operational amplifiers etc. 
{ 5 }. 
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In fact, reference {4} was the only paper on CCD- recursive 
filters presented at the recent conference on CCD-appli- 
cations in San Diego. 

This thesis will study the second order recursive comb filter 
and employ a simulation scheme to study its application 
as a canceller in an MTI-radar. 
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II. SAMPLED ANALOG SECOND ORDER RECURSIVE COMB FILTERS 



A. THEORY 

1 . Analogy Filter 

Since the sampled analog recursive filters are 
implemented by the same circuit configuration used for digital 
recursive filters, we may apply the existing well-developed 
theory for such filters. However, proper modifications must 
be made to account for the fact that the signal is nov; sampled 
analog and that there now are N delay stages instead of one. 

It is known that two general design approaches have been 
developed for the HR (infinite impulse response) type digital 
recursive filter, namely the indirect and direct methods. 

The indirect method starts with the transfer function in the 
LaPlace transform variable s (H(s)), of a continuous 
filter and follov;s by per'forming another transform to obtain 
the corresponding transfer function in the discrete time 
variable z . Several different transforms are available 
for this purpose, for example mapping by differential 
transform, impulse invariant (standard z ) transform, 
bilinear z-transform and matched z-transform. 

The direct m.ethod is digital design in frequency or tim.e 
domain using some type of computer aided design procedure. 

Here the indirect approach will be used: 

Given a prescribed response, there are several methods by 
V7hich one can develop a transfer function in the variable s. 
Examples are: Butterworth, Chebyshev and elliptical filters 

etc{l}. In this thesis the Butterworth approach is used to 
obtain the transfer function in the s-domain. 
a. For a Low Pass Filter 

For the Butterworth filter in the second order 
case the equation is given as: 



where w is the 3-db cutoff frequency. 
Converting to s-domain and expanding: 



H(s)‘ H(-s) 



1 



2N 



Let N =2 



1 + s 




X 



By evaluating the poles in the s-plane and discarding the 
two poles in the right hand plane, one gets: 



2 • Transformation to Z-domain 
a. Low Pass Filter 

By indirect design method, using the bilinear 

z- transform: 

s = 2(z-l) 

Tl?flT 

one gets the following transfer function in z: 



By employing a little algebra one will end up with this 
equation : 



H(s) 



1 



2 

s + 2w s + w 



X X 



X 




H(s) 



2 



s 




X 



2 



2 



w 



H(z) 



X 
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H(z) 



K ( z~~ + a 2 z~^) 

1 + z~^ + b2 z~^ 



where : 



1 



K = 



a 



o 



a-i = 





rp2 2 

T (0 

X 




4 + 2/2 CO T + CO ^ 

X X 


CM 


1 

2 

1 

2 2 

2 CO T - 8 

X 




4 + 2/7 CO T + 


2 „2 
CO T 


X 


X 


4 + CO ^ - 2v^ 


CO T 


X 


X 



4 + 2 /Tw T + CO ^ 

X X 
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b. High Pass Filter 

Starting from the s-domain equation: 



H(s) = 

s + /2 u) s + (0 

X X 



and using the bilinear transform again, one end up vjith: 

-1 - 9 ) 

K( a + a, z + a« z 
H(z) = 2 i 2 

-1 -2 

1 + z + z 



where : 



K 



4 + 2/7 






2 2 



a 



o 



= 1 



a 



1 



-2 



b 



1 



2 2 

2 to T - 8 

X 

4 + 2 /7 to T + w 

X X 



4 - 2/2wT + to 

, _ X X 

4 + 2 /7 to + to 

X X 

It is noted that to^ and T always occur together in the final 
expressions and that there are only two differences between 
the high pass and, lov? pass cases: The size of K, and the sign 
of a^^ . But K is merely a scaling factor, and has no influence 
upon the shape of the frequency response. 

This means that when the coefficients are kept constant , the 
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Coefficient values 



TABLE I 



w T 

X 



w T 

X 



b 



2 



0 






t : — 



.1 
. 5 
.7 



1 . 

1.3 

1.6 

1.9 

2 . 

3 . 

4 . 

6 . 



- 1.86 .5 

- 1.32 1 . 

- 1.09 1.5 

-.77 2 . 

-.49 3 . 

-.26 4 . 

-.06 5 . 

0 6 . 

.47 7 . 

.77 
1.24 



.5 

.28 
.18 
.17 
. 21 
.28 
.34 
.4 
.46 




19 



cutoff frequency will change by varying T. 

Table I and Figure 2.1 shows how b, and b„ varies with w T. 

12 X 

Table II is a sunimary of how the frequency response behave 
as co^T is varied, giving the ranges where the filter acts as 
integrator or canceller. 

Since region I (Figure 2.1) is the area of interest in this 
thesis, the experimental study was carried out here. 

3 . Direct Design Method in Z-Domain 

a. By Trial and Error Method 

By carrying out the design in the z-domain, 
transformations from H(s) to H(z) can be avoided: 




FIGURE 2.2. Z-Domain Unit Circle 

Figure 2.2 shows the z-plane, where the entire left half of 
the s-plane has been transformed into the inside of the unit 
circle. By moving the two poles and two zeros (second order) 
around and using the geometric method {1} to calculate the 
amplitude- frequency response, it is possible to find a 
combination of the singularities that gives a response quite 
close to the- desired response. Appendix A shows the result 
of such a study. 

b. By Computer Approach 

Several computer programs {1} have been 
developed that v;ill optimize the search of coefficients, 
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TABLE II 



Relationship between w T and Filtertype for Highpass 
and lowpass filters using the bilinear z-transform method. 



Filter Type OJ T Filter Shape 




^s = sampling frequency N= number of delay stages 
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based on different approaches in reducing the deviation 
between presented response and filter response. 

This approach was not employed here. 

B. EXPERIMENT 

1 . System Configuration 




FIGURE 2.3. Second Order Recursive Filter 

Figure 2.3 shows the block diagram of the canonical form of 
a second order recursive filter. Serial analog delay (SAD) 
from Reticon Co were used as delay lines. 20 K potentio- 
meters v;ere used to implement the coefficients. In the cases 
where coefficients greater than one were needed for the "a”- 
coef f icients 5 the numerator in the transfer function H(z) was 
scaled down so that the "a"-coeff icients were normalized. 
Several methods to implement the coefficients using the 
potentiometers are discussed in references {2} and {3}. 

In this work Figure 2.4 is presented, giving the calibration 
of potentiometer setting for desired coefficient value. This 
curve was obtained by simply measuring input and output 
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FIGURE 2.4. Potentionieter Calibrat 



across the potentiometer. 

2 . Results 

The clock frequency for the SAD-100 (Reticon Co) 
device was 72 KHZ. Since the effective clock frequency was 
four times 72 KHZ, and there was 96 stages within the device, 
the corresponding nulls of a high pass comb filter should 
occur at multiples of: 

= 3 KHZ 

A computer program (APPENDIX B) was used to calculate and 
plot the theoretical frequency response curve, and also to 
superimpose the actual data points on the same plot. 

The "a"- coefficients were kept fixed, and several combina- 
tions of b^ and b 2 were used. They were selected from 
Figure 2.1, and b^ was negative in all cases studied in this 
thesis . 

By using the values of b^ and b 2 corresponding to w^T, the 
resulting response should be flat-topped in every case. 

The implemented coefficients were: 
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The coefficients and were all picked in Region I 
(Figure 2.1). 

Figures 2.9 through 2.18 show the series of results. 

There seemed to be the same problem with the second order 
filter as was experienced with the first order filter{2}, 
namely a skewness at higher frequencies. Also the results 
show an undesired dip in the passbands, indicating that the 
coefficients are not correctly implemented. 

Suspecting that the experimental implementation inherently 
makes the system frequency dependent, several checks V'/ere 
made for its investigation. The gain of individual delay lines 
was checked for frequency dependence, using both the oscillo- 
scope and a narrowband voltmeter. 

The individual potentiometers input-output characteristics 
were measured when the filter loop was closed. 

It was found that the gains of the delay lines were slightly 
frequency dependent as measured on the oscilloscope, and 
severely and inconsistently frequency dependent when 
measured on the narrowband voltmeter, indicating loading 
effects in the latter case; (the input impedance of the 
voltmeter was 100 Kohm) . 

Figures 2.5 through 2.8 show the calibration of the 
potentiometer settings including the frequency dependence 
for the coefficients used in the previous measurements. 

Figures 2.19 through 2.21 show the fr'equency response when 
the frequency compensated coefficients are used. 

The potentiometer settings v;ere changed every 3 KHZ. 

During several of the "runs" it can be seen that there is 
a slight shift in the position of the experimental curve 
relative to the theoretical curve. Th5.s is due to 
uncontrollable drift of the basic clock frequency in the 
course of the datataking. 

The very latest set of data (Figure 2.21), was taken while 
adjusting a^ and ^2 at every null, after first having set all 
coefficients according to the curves given in Figures 2.5 - 
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2.8, at that null. 

Table III shows the potentiometer ■ settings for a^= 0.5, 
a^= -1.0, a-2= 0.5, b^= -0.7, b2= 0.27: 



TABLE III 





From Curves 


After 


readjus 


ting 


Frequency 


a 

o 


^1 


^2 


^1 


^2 




a 

o 


^2 


3 KHZ 


.714 


-1.0 


. 846 


-.873 


. 407 




.652 


.885 


6 " 


. 714 


-1.0 


.846 


-.873 


.405 




.652 


. 885 : 


9 " 


.714 


-1.0 


.847 


-'.875 


.403 




.689 


.863 


12 " 


. 715 


-1.0 


. 848 


-.876 


.4 




.685 


.859 


15 " 


.718 


-1.0 


.85 


-0.877 


.398 




.685 


.859 


18 " 


.722 


-1.0 


.852 


-.877 


. 396 




. 68 


. 862 



Figures 2.19 through 2.21 do not show any sign of the 
skewness of previous data sets, but the shape of the 
observed curves still deviate significantly from the 
theoretical curves, indicating that the proper coefficients 
have not been found for all frequencies. 
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FIGURE 2.5. Frequency Variation of a =.5. a^=.5 and 
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FIGURE 2.6. Frequency Variation of = .3, and 
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FIGURE 2.7. Frequency Variation of b 
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FIGURE 2.8. Frequency Variation of b, = .2, b«= .19 and b«= .21. 




FIGUR2 2.9a. Theoretical and Observed Freqxiency Response of High Pass 
Filter; = -]..0 , = 0.38 
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FIGUaS 2.9b. 



Theoretic^^l and Observed Response of High Pass ?’ilter; 



= -1.0 , b2 = 0.38 
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Magnitude (nonraliaed) 



figure 2.10a. Theoretical and Observed Response of Loi^ Pass h liter 
=■ 1.0 , = 0.38 
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Magnitude (normlized) 



FIGUdE 2. 10b. Theoretical and Observed Response of Low Pass Filter; 

- 0.38 
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FIGUre 2.11a 



Theorstical and Observed jccsponse of High iaas Filter; 



= -0.7, = 0.27 
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FIGUlJi 2, 11b. Theoretical and Observed Response of High Pass Filter; 

= -0.7, bg = 0.27 
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llagnitude (normalized) 



FIGUPJ5 2,12b. Theoretical and Observed Respon-se of Low Pass Fj_lter; 
bi = 0.7, b2 - 0.27 
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Theoretical and Observed Response of High Pass 
- -0.4, “ 0.21 



Filter; 



39 




FIGURS; 



2.13b. 



Theoretical and Observed Response of High Pass Filter 

b, = -o.U, b^ = 0.21 
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FIGUiffi 2.14a. 



Theorotioal and Observed Renponss of Low Pass PUter 

= 0 . 4 , = 0,21 
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FIGUHJ.'^ 2.14b. Theoretical and Observed Response of Lov/ Pass Filter; 

= 0.4, bg= 0.21 
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FIGUFlE 2.15a. 



Theoretical and Observed Response o.f High Pass Fxlter 

= -0.3, ^2 = 0.2 
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FIGUfiJ5 2,15b. Theoretical and Observed Response of High Pass Filter 

bg - 0.2 



= - 0.3, 




FIGURID 2.16a, Theoretica.1 and Obs(5rved Response of Lov-f Pass Filter 

= 0,3, = 0.2 
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Theoretical and Observed Response of High Paso Filter, 
- -0.2, b2=0.19 
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FIGURK 2,17b. Theoretical and Observed Response of Hi"h Pass Filter; 

= - 0.2, b^ = 0.19 
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FIGURii 2.18a. 



Theoretical and Observed .iesponse of W./ Pass Filter; 

b ^ “■ 0,2^ bg ~ 0 , 1 9 
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FIGU.’u:: 2.18b. Theoretical and Observed Response of Low Pass Filter; 

^ 0.2, = 0.19 
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FIGURF. 2.19a. 



Theoretical and Observed Flesponse of Hi^h pass filter 

= - 1 . 0 ; = 0-38 
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FIGUi'u^ 2.19b. Theoretical and Observed Response of High Pass Filter; 

= “1 . 0 . b ^ = 0.38 
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FIGUitJC 2.20a. 



Theoretical and Observed Response of I.ov Pass Filter 

^^-+1.0, b2 = 0.30 
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FIGU'U*: 2.20b. Theoretical and Observed Response of Lov/ Pass Filter 

= -UO, - 0.38 
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= - 0 . 7 > ^2 ~ 0,27 



III. APPLICATION TO MTI-RADAR 



A. THEORY 

1 . Statement of the Radar Signal Processing Requirement 

The majority of existing radars are pulsed radar 
systems. Thus the radar return is in effect a sampled signal, 
and the sampling width is the pulse width (PW) of the radar. 
Figure 3.1a. is a tim.e-amplitude plot of radar signals. 

At multiples of T= the transmitted pulses are shown. 

The time interval between each transmitted pulse, called 
the Pulse Repetition Interval (PRI), target returns will 
appear. In Figure 3.1a two stationary targets and one moving 
target are shown. The stationary target return may be either 
positive or negative depending upon the phase of the returned 
signal. There is, however, no change from interval to interval, 
a fact whicli is exploited in the MTI canceller. 

The moving target however, is changing in amplitude every 
interval, at a rate which is the doppler frequency. 

Figure 3.1b shows how the pulse train of Figure 3.1a 
would look as presented on an oscilloscope or A-scope. 

Here every PRI is shown on top of each other, and the moving 
target is exibiting the socalled "Butterfly effect" {6} 
while the stationary targets remain fixed in amplitude and 
time . 

To get the frequency representation of a periodic pulse 
train in the time domain, one can take the Fourier transform.. 
The result v;ould be a sin x -type waveshape. If only one 

X 

pulse were transmitted, the fr'equency spectrum would be 

continuous, and the first null of the curve would appear at 
1 

“ PW* But wlaen transmitting a train of pulses spaced in time 
at T - , ideally each target will return a square pulse 

at PRI equal to T. frequency sx->ectrum of this periodic 

signal can be found from the Fourier series if the pulse 
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FIGURE 3.1. Infinite periodic Pulsetrain (a), and its appearance on an 
Oscilloscope (b). 



train is infinite: 



jnoj- 1 

f(t) = Z C e ^ 
n 



where 



C 



Y~ f(t) e 



T/2 -jnw^t 



dt = 2n 

T 



n 



2 



is shown in Figure 3.2, where the sin x modulation 

X 

whose zero crossings with co- axis are determined by the 
pulsewidth. 

Example: For a radar with PRF of 800 HZ, a PW of l.Sysec 

and v;avelenght X of 0.1 meter. 

The first zero crossing of the envelope will occur at 
769 KHZ. For a target radial speed (v^) of 900 knots, 
using a formula from ref {6}; Doppler frequency, f^= 2v^ 



Thus it is obvious that the peaks of the repetitive 
"clutter spectrum" are essentially of equal amplitude in the 
frequency range of interest in MTI-radar (Figure 3.3). 

Figure 3.3 is a small portion of Figure 3.2 (near the origin) 
and one moving target is introduced in this ideal frequency 
spectrum which consists of impulses. 



X 




2 • 900 . 1852/3600 

0.1 



9.26 KHZ 
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FIGURE 3.2. envelope of many impulses at 

f = n • PRF 



Figure 3.4 shows the "real" clutterspectrum, i.e. the 
impulses have been spread out in frequency. There are 
several reasons for the nonideal shape shown; The energy 
in the returned echoes is spread out in frequence due to 
land target movement (vegetation swaying in the wind), sea 
waves moving, platform noise (own irregular movement), 
and effects of a finite number of radar returns which stems 
from the fact that a radar antenna normally rotates. 

Thus v;e do not have an infinite sequence of target pulses. 
Rather, the number of "hits" on a particular target may be, 
typically, from 10 to 50 per antenna rotation period. 

Figure 3.5 shov;s the bursts of target echoes every time the 
antenna points at the target; the time between each burst 
is called the antenna rotation period. 

In Figures 3.6 and 3.7 the Fourier transform of a finite 
number of 6, respectively 10 hits are shown. The mathematics 
is omitted here. 

The fewer the number of hits, the more energy leaks into 
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FIGURE 3.3. Ideal Frequency Spectrum of Clutter and One Moving Targe 




FIGURE 3.4. A "Real" Clutter Spectrum 



frequencies where it may give raise to false alarm. 

Thus the finite number of hits not only leads to spreading 
of the clutter spectrum around multiples of PRF , but also 
reduces the maximum attainable cancellation ratio. For 
proper MTI operation, at least 10 hits per scan is needed. 

The Fourier transform of such a truncated series of target 
returns in which there is a moving target among heavy clutter 
return, is shown in Figure 3.8. Note that the component is 
also spread out from the ideal impulses, and the energy in 
the sidebands. 

b. Tim.e Response Approach to MTI-Radar Analysis 

In this analysis we do not concern ourselves v?ith 
the trailing edge of a typical echopulse, as our filter has 
no noticeable disadvantageous effect in this area. 

But since our filter is a recursive filter, the resirculation 
of signals lead to undesirable effects v;hich must be 
investigated , 
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FIGURE 3.5. Finite Bursts of Target Echoes 
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FIGURE 3.6. Fourier Transform of Periodic Bursts of 6 Pul 
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FIGURE 3.7. Fourier Transform of Periodic Bursts of 10 Pulses. 



4-1 
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FIGURE 3.8. One Moving Target Among Heavy Clutt 



MT I -video 

>• 




FIGURE 3.9. Delay Line Canceller 



Figui'^e 3.9 shows the blockdiagram of a conventional MTI 
delay-line canceller {G}. Cancellation is obtained only v:hen 
two succeeding pulses are equal in amplitude. The first and 
last pulse in a pulse train of a finite number of pulses will 
be passed by this filter. 

In the recursive CCD-filter of Figure 2.3, the first 
pulse of such a fixed amplitude pulse train will be passed, 
the second pulse will not be fully cancelled. After the 
pulse train input is terminated, it will take a finite time, 
depending upon the values of the feedback coefficients b^ 
and b2, before the filter output approaches zero. 

This time is measured in number of PRI's, and the effects in 
PPI presentation is that the target will spread out in 
azimuth in the direction of the antenna rotation. 

Figures 3.10 and 3.11 illustrates typical unit impulse 
responses of a first order recursive filter v;ith feedback 
coefficients equal to 0.3 and 0.9, respectively. 

Example: A typical antenna rotation rate is 20 revolutions 
per minute, which means 3 seconds/scan or 8.33 msec/degree. 
For the case where b^ = 0.9 this means that the target 
would be extended 29/8.33 = 3.5 degrees beyond the actual 
position. This is considered to be a serious limitation, 
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Anpl (norm) 




FIGURE 3.10. Time Response of First Order Recursive 
Filter^ = 0.3 




FIGURE 3.11. Time Response of First Order Recursive 
Filter, = 0.9 
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and consequently we have to keep this aspect in mind when 
designing filters of this type. 

2 . Applying the Comb Filter Technique. 

The purpose of MTI -radar is to detect targets 
moving at or above a certain minimum speed, and suppress all 
echoes from stationary targets (clutter). 

To do this, one exploits the well known doppler effect, 
namely that the frequency of a returned echo from such a 
target (moving), is different from the transmitted frequency. 
Using def ectioh~"techniques in the~recoiver {6} , {7}, only 
doppler frequencies are“amplif ied and become the KTI- 

video which may be transmitted to a PPI for viewing. 

Thu^ we need a_comb filter (Figure 3.12), which can reject 
the clutter, and pass the dopp-ler shifted signals. 

One disadvantage is that when the target's I'-adial velocity 
corresponds to some multiple of the PRF, it will be 
suppressed also (the target is then said to have a radial 
blind speed). To overcome this one must var’y the PRI fi^om. 
pulse to pulse, by PRF-jitter (random variation of PRI), 
or PRF stagger (sequentially stepping the PRI between two 
val.ues only). 

A proposed method for building a CCD - filter to facilitate 
such-a system is outlined in chapter IV. 




PRF 2 PRF f 



FIGURE 3.12. Comb Filter and Clutter 
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3 . Bandwidth Requirements. 



The pulsewidth determines the bandwidth (BW) 
necessary to reproduce the signal. 



BW 



1 

= 2W 



The employed delay lines have 96 stages each. The necessary 
clock frequency, f^, is then: 



f = PRF- 96 
c 

For PFxF = 800 HZ, ~ 7 6.8 KHZ, which is the effective 

sample frequency (f^). 

Using Nyquist's Sampling Theorem, f = 2BW. 

s 

The maximum, bandv/idth that can be faithfully reproduced is 
then 38.4 KHZ. This corresponds to PW = 13 ysec, which is 
quite v;ide and unrealistic. Thus delay lines with considerable 
larger number of stages are needed. 

B. EXPERIMENT 

1 . Simulation of Radar Return Signal 

Since the size of the delay line v/as limited to 9 6 
stages it was decided to simulate the radar signal for which 
the delay line was able to respond. The objective was to 
provide radar return signals simulating fixed and moving 
tar’gets with variable PRF and number of hits per scan. 

The simulation of PRF was implemented using a Beckman 9030 
function generator controlled by a sinewave generator 
(Wavetek Model 142). 

The function generator produces a square pulse whose lenght 
and amplitude can be varied for every period of the sinewave. 
Thus the PRF may be controlled by changing the VJavetek 
frequency (Figure 3.13). 
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"FIGURE 3.13. Time Relationship Between Output from 

Function Generator and Sinewave Generator 



The pulse train thus obtained is the simulation of clutter, 
or fixed targets, since each pulse is equal in amplitude. 




'Wavetek 



FIGURE 3.14. Block Diagram of Simulation Configuration 




Figure 3.14 illustrates the simulation setup. A National 
Semiconductor MM 4616A switch (Figure 3.15) was used to 
perform the switching functions. The specifications of this 
arrangement is that each pulse through this switch had a 
different amplitude, simulating the moving target signal. 

On the oscilloscope, this signal should exhibit the "butter- 
fly effect" {6}. If needed, both stationary and moving target 
signals could be fed into the filter simultaneously. 

The second switch was operated by another Wavetek M142, which 
was set such that the positive part of the squarewave was 
much shorter than the negative part. The positive part 
would turn the switch on, letting a finite number of pulses 
from switch one through. This was the simulation of the 
antenna sweep. However, the present simulation of this 
effect is limited because it has not been possible to 
simulate the long time between two series of radar hits. 

2. Results. 



a. Stationary Targets with PRF = 800 HZ 




FIGURE 3.17. Cancelling of "Stationary" Target, 

Horizontal: 10 ysec/cm. Vert: .5V/cm 
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FIGURE 3.15. Analog Switch Schematic. 
TABLE IV 

Typical Radar Charachteristics 



Band: 


j-band 


S-band 


C-band 


X-band K 


-band 1 


Frequency 

(GHZ) 


L - 2 


2-4 


4-8 


8-12.5 1 


2.5-18 


Typical 

PRF 


200- 

500 


300- 

600 


500- 

1000 


800- 2 
3000 


000- 

5000 


PRI= 1 
PRF 


2-5 

msec 


1. 67- 
3 msec 


1-2 

msec 


333- 

1250 usee 


200- 

500 usee 


Maximum 
Theoretical 
Range n.m. 


L62- 

405 


135- 

270 


81- 

162 


27- 

101 


16.2- 

40.5 


Typical PW 


3-5 

usee 


1-3 

usee 


0.5- 
4 usee 


.05- 

1 . 5 usee 


.01- 

0.5 usee 
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FIGURE 3.18. Cancelling of Stationary Target. 

PRF= 800HZ. Horizontal; 10 usec/cm. 




FIGURE 3.19. Cancelling of Narrow, Stationary Target. 

PRF = 800HZ. Horizontal: 10 usec/cm. 
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Figure 3.20 shows a "moving target" signal. 'The doppler freq- 
uency was at the center of first comb tooth of the filter. 
It can be seen that the combfilter did not cancel the signal. 
The insertion loss was quite small. 




FIGURE 3.21. Moving and Stationary Targets. PRF= 800 KZ 
Horizontal: 10 ysec/cm 
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c. Both Stationary and Moving Signals, PRF = 800 HZ 

Figure 3.21 shows an input consisting of both moving 
and stationary target signals. It can be seen that only the 
moving signal is passed. 

d. Fixed and Moving Targets. PRF = 10 KHZ 

To try to eliminate some of the glitches, the effect- 
ive clock frequency was increased to 288 KHZ, so that the 
nulls of the comb filter appeared at multiples of 10 KHZ. 

The set of coefficients for the following results were 
(potentiometer settings): 

a = 0.6885 , a^ = 0.975, a„ = 0.87, b, = 0.5475 , b.,= 0.498 

O 1 2 ’1 2 

For all the pictures the vertical scale is 0.5 V/cm. 




FIGURE 3.22. Stationary Target Signal. PRF = 10 KHZ. 
Horizontal scale: 50 ysec/cm. 



Figure 3.22 corresponds to Figure 3.18, the difference 
being that now the PRF was 10 KHZ instead of 800 HZ. Note 
that the glitches were much smaller in amplitude at higher 
clock frequency. 
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FIGURE 3.23. Moving Target Signal only. PRF = 10 KHZ 
Horizontal scale: 50 ysec/cm. 



e. Moving Target. PRF = 10 KHZ. 

Figure 3. 23 shows a moving target only, with 
doppler frequency, = 5 KHZ. Again it is noted that the 
moving target signal passed through the comb filter with low 
insertion loss. 




FIGURE 3.24. Moving Target, f^ = 8.5 KHZ. PRF+ 10 KHZ. 
Horizontal: 50 psec/cm/ 
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f. Moving Target with Doppler Frequency Close to the 
First Null 

Figure 3.24 shows that the output amplitude indeed 
was attenuated when was close to the blind speed. 

In this case the first blind speed was the radial velocity 
of the target corresponding to f^ = 10 KHZ. 




FIGURE 3.25. Moving Target and Stationary Target, 
f^ = 5 KHZ, PRF = 10 KHZ. 

Horizontal: 50 ysec/cm. 

g. Stationary and Moving Targets, PRF = 10 KHZ 

I 

Figure 3.25 shows that when the input consisted of 
both stationary and moving targets, only the moving target 
signal was passed. In this case f^ = 5 KHZ which was midband 

h. Stationary and Moving Targets, Doppler Frequency at 
First Null 

Figure 3.26 depicts the same situation as in Figure 
3.25, only now f^= 10 KHZ which was the first blind speed. 
Both targets were suppressed, and the glitches are seen. 

The output was almost identical to the output in the station 
ary target case (Figure 3.22) 
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FIGURE 3.26. Moving and Stationary Targets Superimposed. 

= 10 KHZ, PRF = 10 KHZ. Horizontal: 50 
ysec/cm. 



i. The Effect of a Finite Number of Stationary Target 
Signals . 




FIGURE 3.27. Truncated Series of Fixed Targets (15 hits) 
Horizontal: 0.5 m.sec/cm. 

Figure 3.27 shows the input and output of the comb filter in 
the case where Switch 2 (Figure 3.14 ) was operated by a low 



78 



frequency asymmetrical squarewave , to simulate finite number 
of "hits" on the target during one period of the antenna scan 
The input to the switch was an infinite train of fixed 
cimplitude pulses. The oscilloscope was synchronized to the 
switch controlling waveform. Note that the first and last 
pulses in each truncated series of radar hits were passed by 
the comb filter. This was in agreement with the time analysis 
The rest of the pulses were cancelled, except for the 
glitches, which appeared quite pronounced in this picture. 
Additional study is necessary to determine the cause of the 
glitches and their elimination. 




t 

FIGURE 3.28. One Pulse from Figure 3.27. 

Horizontal scale; 10 ysec/cm. 



Figure 3.28 shows one of the individual pulses in the 
truncated trains of radar hits from Figure 3.27. The 
oscilloscope was now synchronized to the function generator 
(Figure 3.14). It is noted that the cancellation was poor, 
verifying the analysis of section III. A. 1. a, which predicted 
that energy would leak into sidebands when a finite number 
of hits were being processed, causing false alarm. 
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FIGURE 3.29. Truncated Series of Fixed Targets (3 Hits). 

Horizontal: 20 ysec/cm. 

Figure 3.29 shows a truncated series of only 3 hits. The 
oscilloscope was now synchronized to the squarewave generator 
controlling Switch 2 , and the three hits cannot be seen 
clearly on this picture for that reason. Note that the first 
two pulses were passed and the last pulse was also passed as 
well as several recirculated pulses, causing a smearing of 
the trailing edge. This phenomenon was desribed in section 
III.A.l.b., and caused by recirculation of the signals. 

The signal will die out gradually depending upon the size of 
the feedback coefficient. 

Figure 3.30 shows one of the individual pulses shown in 
Figure 3.29. Although this pulse was similar to Figure 3.28, 
it should be quite different. More pronounced false alarm 
was expected in Figure 3.30 than in Figure 3.28. The reason 
for the poor result here was that the individual truncated 
series were too close, so that several consecutive series 
were superimposed on the scope. This was also the reason 
for the base line to appear even under the pulse in both 
Figures 3.28 and 3.30. 
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FIGURE 3.30. One pulse from Figure 3.29. 

Horizontal: 10 ysec/cm. 

The relative brightness of the pulse itself and the base-line 
under the pulse in Figure 3.28 indicated that there were more 
oscilloscope sweeps for pulses than for base line. 

But in Figure 3.30 the situation was reversed. This v;as 
reasonable, since here most of the time-axis was a base-line. 






81 




A. CONCLUSIONS OF EXPERIMENTS 



1 . Filter Experiments 

It has been shown that the CTD-device can be used 
for implementing second order comb filters. 

2 . MTI-Radar Simulation Experiment. 

Most of the simulation objectives were met except for 
the simulation of antenna rotation. Due to the limitation of 
instrumentation, the time between two series of radar hits 
is too short. A more realistic case would be to have a series 
of radar hits appear every two seconds. 

A major unresolved problem encountered in this study was the 
observed glitches. It must be solved before the sampled 
analog recursive CTD comb filter may be successfully employed 
in real radar systems. 

Characteristics of the present CTD-devices do not yet allow 
for enough bandwidth. More delay stages are needed. 

Until then only simulated results may be obtained, with 
rather restrictive ranges of pulsewidths and PRF’s. 

However, deviations between theory and measured results have 
been noticed. It was found that both the delay line and its 
peripheral electronics do not have adequate bandwidths . 
Consequently, the coefficients used in experimental studies 
vary V7ith frequency. It is believed that this is one of the 
major causes of deviation of experiment from theory. 

B. RECOMMENDATIONS 

1 . Research Topics Directly Related to Present Thesis. 

In the filter study of chapter II, more work is 
needed to achieve implementation of frequency independent 
coefficients. Perhaps voltage controlled operational 
amplifiers should be used instead of potentiometers, and 
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also a frequency independent delay-line and its peripheral 
circuitry. 

In the MTI- simulation study, the problem is the glitches. 

To find the causes, and to eliminate them is important. 

Another area is to build a better simulation setup for 
examining the effect of finite number of returns. 

2 . Proposed Application to Staggered/ Jittered PRF MTI- 
Radar . 

Not too many military MTI-radars today will be 
built that do not include either jitter PRF or- stagger PRF 
capabilities. One of the major attractions of CTD-delay 
lines is the ease and speed with which the clock-frequency, 
and thereby the delay time, may be changed. Consequently, 
the null frequencies of the comb filter can be easily changed 
by changing the clock frequency. 

But there are major problems in applying the CTD-filter to 
such cases. This is because in jitter, or staggered PRF -radar, 
the PRI changes from pulse to pulse, and thus the clock 
fi’equency of the delay line must be changed for every PRI. 

•It is seen that information which has already been fed into 
the delay-line (memory device), will be clocked out at a 
different clock frequency and then added in the output summer 
(Figure 2.3) with a signa employing a third PRI. The result 
would be meaningless. 

Figure U.l is a proposal to employ clocked CTD-delay lines 
in a jitter PRF tracking filter which could overcome this 
problem. It may solve the undesirable effects of long 
settling time associated with high feedback coefficients, 
as stated in section III.A.l.b. 

A short explanation using difference equations, of this 
circuit follows: 
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Counter reset 
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FIGURE 4.1. First Order CCD Jitter ?RF MTI-Filte 





1st pulse arrives: 



Upper part of 
filter effective: 

Video into lower 
circuit, still 
tapping upper part: 



SI opens. Video into 
upper part : 

SI closes, S2 opens. 

Video into lower part: 

SI opens, S2 closes 
Video into upper part: 

I 

If the video train input is 
pulses, then: 



Y(0) = a^ X(0) 

Y(l) = a^X(l) + (bjLa^+a^)X(0) 

Y(2)=a X(2)+(bTa^+aT )X( 1 ) + 
o 1 o 1 

bj^^a^X(O) + a^b^X(O) 
Y(3)=a^X(3) + (b^a^+ a^)X(2) 

Y(4)=a X(4) + (b,a + a^ ) X(3) 
o 1 o 1 

Y(5) = a X(5) + (b^a +a, )X(4) 
o 1 o 1 

etc 

a train of fixed amplitude 



Y(n) + X(n){(a (l+bT))+aT} 
oil 

There will be perfect cancellation if: 

a^(l + b, ) + ai = • 0 
o 11 

Circuit operation: 
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Vin ^ sawtooth waveform. A sample and hold circuit is used 
to establish a DC-voltage level which determines the 
frequency output of the voltage controlled frequency genera- 
tors. Starting with frequency generator if 1, assuming that 
the video is being fed into the upper part of the filter, the 
binary counter will count 96 periods and then generate a 
trigger pulse, which controls the radar transmitter and all 
radar timing circuits and the various switches in the filter 
itself. The functions of this trigger pulse in the filter are 
to: 

- Reset the lower counter to zero. 

- Sample V. to establish a new f « . 

- Switch the video to the lower j^iart of the filter. 

- Switch off S2 , the output of the lov7er delay line. 

- Switch on SI, the output of the upper delay line. 

Next, after another PRI and f ^2 counted to 96, a trigger pulse 
appears (the counters will only produce a trigger pulse on the 
count of 96 after each reset), f^^ is then set to a new value, 
counter ^ 1 is reset, the video is switched to the upper part 
of the filter, the switch SI is opened and S2 is closed. 

Thus the direct signal will always be summed in one part of 
the filter with the delayed signal stored in the other part. 
This scheme may be extended to higher order filters. 
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APPENDIX A 



FILTER DESIGN IN THE Z - DOMAIN 

It was desired to design a filter with 100 HZ across 
the null at an attenuation of at least 20 db down. 

From a pole-zero plot in the z-domain Figure A-1), the 
magnitude of the transfer function can be evaluated for every 
point on the unit circle, geometrically. 

By changing the pole and zero locations within the unit 
circle in a trial and error manner, and comparing the 
magnitude of the transfer function with the desired shape of 
the frequency response curve, it was possible to come fairly 
close after a rather short time: 




I 



FIGURE A-1. Pole-Zero Plot in Z-Domain. 



Results: Zeros: Z ^2 “ - j0*l 

Poles: P ^2 = 



This gives the transfer function: 



^ ' O.SSG - z“^ + 0.4556 z" 
H(z)= ^ 

. 1 + 0.14 z" + 0.2074 z 
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A computer simulation of' a filter using these cofficients 
is shown in Figure A-2. The desired width of the stopband, 
and attenuation was obtained. The filter was also built, 
and a few data points were plotted, but the implementation 
of the coefficients were not correct. 
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FIGURE A-2. Frequency Response of Filter Using Z-Domain Design 
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